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ABSTRACT 


Consider  a reliability  system  that  is  composed  of  n components 


each  of  which  is  operating  at  some  performance  level.  We  suppose  that  there 

exists  a nondecreasing  function  4>  , called  the  structure  function,  such 

. '■  > 

that  >)>(x^,  x^)  denotes  the  performance  level  of  the  system  when  the 

• th  ^ ' , , . . , 

1 — component  s performance  level  is  x^  , i = 1 n . ..  t 

Whereas  almost  all  previous  work  assumed  that  both  x^  and 
(}i(x^,  x^)  were  binary  variables  we  shall  allow  both  to  be  arbitrary 

nonnegative  numbers  and  we  extend  many  of  the  important  results  of  the  usual 
binary  model  to  this  more  general  framework.  In  particular  we  obtra-i-B  a 
fundamental  inequality  for  E[<()(X^,  ...,  X )]  when  (()  is  binary  which  can, 

r' 

among  other  things,  be  used  to  generate  a host  of  inequalities  concerning 
IFRA  distributions  including,  as  a special  case,  the  IFRA  convolution  theorem. 
In  Section  2 we  define  the  concept  of  an  IFRA  stochastic  process  and  prove 
the  analog  of  the  IFRA  closure  theorem;  and  in  Section  3 we  do  the  same  for 
NBU  stochastic  processes.  In  the  final  section  we  present  some  applications 
to  stochastic  networks. 
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MULTI-VALUED  STATE  COMPONENT  RELIABILITY  SYSTEMS 


by 

Sheldon  M.  Ross 


0.  INTRODUCTION  AND  SUMMARY 

Consider  a reliability  system  that  is  composed  of  n components 
each  of  which  is  operating  at  some  performance  level.  We  suppose  that 
there  exists  a nondecreasing  function  ij)  , called  the  structure  function, 
such  that  (}>(Xj^,  x^)  denotes  the  performance  level  of  the  system 

when  the  i—  component's  performance  level  is  x^  , i = 1 , . . • , n . 

Whereas  almost  all  previous  work  has  assumed  that  both  x.  and 
(()(x^,  x^)  were  binary  variables  we  shall  allow  both  to  be  arbitrary 

nonnegative  numbers.  In  the  next  few  sections  we  extend  many  of  the  important 
results  of  the  usual  binary  model  to  this  more  general  framework.  In 
particular  we  obtain,  in  Section  1,  a fundamental  inequality  for 
E[(}i(X^,  ...,  ] when  <p  is  binary  which  can,  among  other  things,  be  used 

to  generate  a host  of  inequalities  concerning  IFRA  distributions  including, 
as  a special  case,  the  IFRA  convolution  theorem.  In  Section  2 we  define 
the  concept  of  an  IFRA  stochastic  process  and  prove  the  analog  of  the  IFRA 
closure  theorem;  and  in  Section  3 we  do  the  same  for  NBU  stochastic  processes. 
In  the  final  Section  we  present  some  applications  to  stochastic  networks. 
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THE  STRUCTURL  FUNCTION 

Suppose  now  the  perforunance  level  of  component  i is  a random  variable 
having  distribution  where  F^(x)  = P{X^  > x}  , and  suppose  that 

the  X^  are  independent.  We  define  the  function  r(Fj^,  ....  F^)  by 

r(F^,  ....  F^)  = E[$(X^,  X^)] 

and  call  r the  reliability  function  of  the  system.  It  immediately 
follows  from  the  monotonicity  of  4>  that 

Proposition  1: 

If  F.  and  G.  are  distributions  such  that  F.(x)  > G (x) 

11  1—1 

for  all  X then 


K>  i ■■<=1 • 


We  shall  need  the  following  lemma  which  is  a slight  variation  of  a lemma 
used  by  Block  and  Savits  [3]  to  prove  the  IFRA  convolution  theorem. 

Lemma  1: 

Let  r(s)  be  a nonnegative  nondecreasing  function  of  s , s ^ 0 
and  let  G be  a distribution  function  with  G(0)  = 1 . Then  for  0 < a 1 , 

(r(s))“d(l  - g“(s))  > 

0 


00 


r(s)d(l  - G(s)) 


Proof: 


I 


Lemma  4.1  (p.  217)  of  ’oss  [5]  (also  given  as  Lemma  2.3  (p.  84)  of 


Barlow  and  Proschan  [1]  generalizes  to  give  that  for 


n 


From  this  the  conclusion  follows  from  a standard  limiting  argument 
(as  in 

The  following  theorem  is  of  fundamental  importance. 

Theorem  1; 

If  (()  is  a binary  function  then 


(*) 


. F“)  > [r(F^, 


F^)r 


for  all  0 < a < 1 . 


Proof : 


The  proof  is  by  induction  of  n . When  n = 1 it  follows  from  the 

1 1 X ^ c 

monotonicity  of  $ that  it  must  be  of  the  form  (Ji(x)  = for  some 

c . Hence  E[())(X^)]  = Fj^(c)  and  so  both  sides  of  the  inequality  (*)  are 
equal.  So  assume  (*)  for  all  binary  structures  of  n - 1 components,  and 
consider  the  n component  case.  Conditioning  on  yields 


(1) 


(F“,  ....  Pt-  JrjF» f;;.j  )<i(l  - f“(s)) 


where 


rs(F“,  ....  = E[<*>(Xj^,  ....  X^_j^,s)] 


with  X^  having  distribution  Fj^  . By  the  induction  hypothesis  we  see  that 


A 


Fi, 


K-i)  i 


••• 


and  so  from  (1) , 


r(Fj,  •■■,  5“)  i ■ 

As  it  follows  from  the  monotonicity  of  4»  that  r^  is  nondecreasing  in 
s we  can  apply  Lemina  1 to  the  above  to  obtain  that 

r "1^ 

■■(fj O i I 


= (r(F^,  ....  F^)) 


Definition: 

The  distribution  function  F , with  F(0)  = 1 , is  said  to  be  an 
IFRA  distribution  if 

F(ax)  > f“(x)  for  all  0 a < 1 , x ^ 0 . 


Corollary  1: 


If  X^,  •••>  X are  independent  random  variables,  each  having  an 
IFRA  distribution,  then  for  all  nondecreasing  binary  functions  <t>  , 


E 

X^,  ... 

’ V 

a 

a 

a 

> (E[4>(X,, 

•••  9 

X )])“ 

for 

Proof : 


If  is  IFRA  with  distribution  F^  , then  P > x^  i 
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and  so  from  Proposition  i 


>(X^, 


...  X^) 


...  f: 


and  the  result  follows  from  Theorem  l.| 

The  above  Corollary  provides  a host  of  inequalities  concerning  IFRA 
random  variables.  For  instance  we  have 


Corollary  2; 

If  X ....  X are  independent  IFRA  random  variables  then 
1’  n 


(a)  I X.  is  IFRA. 
1 ^ 


(b)  P jn  X.  > aa'^j  > |pjn  X.  > a|j 


0 < a < 1 . 


Proof : 


Part  1 follows  from  Corollary  1 by  using  the  function 


to  obtain 


jl  if  I X > a 
1’  n^  i 

lo  otherwise 


(pjlx.  >aj) 


similarly  (b)  follows  by  using 


|i  if  nx.  2^ 


’(X.  > • • • » X ) “ j 

^ (O  otherwise.! 


I 
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2.  THE  GENERALIZED  IFRA  CLOSURE  THEOREM 

In  this  section  we  suppose  that  the  component  performance  levels 
vary  with  time  and  we  let  X^(t)  denote  the  level  of  component  i at 
time  t . Thus,  for  instance,  <f(X(t))  = <()(Xj^(t),  ...,  X^(t))  denotes 
the  systems  performance  level  at  time  t . 

Definition; 

The  real-valued  stochastic  process  {X(t),t  ^0}  is  said  to  be  an 

IFRA  process  if  T is  an  IFRA  random  variable  for  every  a , where 
^ a 

T = inf  {t  : X(t)  < a} 
a — 

is  the  first  time  the  process  reaches  or  goes  below  a . 

Theorem  2:  The  Generalized  IFRA  Closure  Theorem 

If  {X^(t)}  ,1=1,  ...,  n , are  nonincreasing  independent  IFRA 
processes  then  {((i(X(t))}  is  also  IFRA  whenever  4>  is  nondecreasing. 


Proof ; 

Let  F (x)  = P{X  (s)  > x}  , and  suppose  first  that  4>  is  a binary 
i , s i 

function.  Let  T denote  the  first  time  t that  (})(X(t))  = 0 . Now 


(2) 


P{T  > at}  = Pt({i(X(at))  = 1}  by  monotonicity 


= E[^>(X(at))] 


F . ) . 
n,at 


I 

I 


Now 
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r.  (b)  = Pix.Cit)  > b} 

i.at  1 

= P{T,  , > at} 

1 ,b 


where  T,  , denotes  the  first  time  that  X.(t)  hits  or  eoes  below 
1 ,b  i 

Hence  from  the  IFRA  hypothesis  on  (X^(t)}  we  see  that 


P{T.  , > at)  > (P{T.,  > t}) 
i,b  — lb 


F.  (b)  > f“  (b) 
i,at  - i,t 


and  so  from  (2)  and  Proposition  1 


PIT  > at)  

> (r  (F^  ^ , . . . , F ))°‘  by  Theorem  1 
~ i , t n , t 


= (P{T  > t}) 


which  proves  the  result  when  (jj  is  binary.  For  an  arbitrary  nondecreasing 
$ we  can  show  that  the  time  to  go  below  b is  IFRA  by  using  the  result 
in  the  binary  case  on  the  binary  function  defined  by 


I 1 if  (()  (}^)  > b 

\ 0 if  4>  (x)  < b 


When  0 is  a binary  function  we  usually  say  that  the  system  fails 
at  the  first  t such  that  <(i(_X(t))  = 0 . Thus  if  each  component  process 
has  the  property  that  the  time  it  takes  to  reach  or  go  below  any  given 


& 


level  is  an  IFRA  random  variable  then  so  is  the  time  to  system  failure. 
The  following  are  two  examples  for  which  the  component  processes  are 
IFR.1. 

Example  1:  A Semi-Markov  IFRA  Process 

If  tX(t),t  2.  ^ ^ is  a semi-Markov  process  such  that  X(0)  = m , 

P.  . , = 1 , i > 1 , P__  = 1 , and  the  time  in  state  i is  an  IFRA 

1,1-1  — 00 

random  variable,  then  it  follows  from  the  IFRA  convolution  theorem 
(Corollary  2)  that  {X(t);  is  an  IFRA  process.  This  would  be  a model 
for  a component  that  gradually  went  to  lower  states  until  it  died 
(reached  state  0),  spending  an  IFRA  amount  of  time  in  each  state. 

Example  2:  A Poisson  Shock  Model 

Suppose  the  component's  level  remained  constant  between  times  of 
extreme  stress  which  occurred  in  accordance  with  a Poisson  process. 

If  at  these  moments  it's  level  decreased  according  to  a given  distribu- 
tion then  the  component  process  would  be  IFRA.  That  is,  if 

N(t) 

C - I X 
i=l 

where  X^  , i ^ 1 are  independent  and  identically  distributed  (i.l.d.) 
nonnegative  random  variables  that  are  also  independent  of  the  Poisson 
process  N(t)  , then  (X(t)}  is  IFRA.  This  result  was  first  proven  by 
Esary,  Marshall  and  Proschan  [4]  who  also  showed  that  the  same  result 
could  be  obtained  under  weaker  conditions  then  the  l.i.d.  assumption  on 

i ■ 


the  X 
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3.  AK  NBU  CLOSURi:  THEOREM 

We  start  with  some  definitions. 


Def inition : 

The  distribution  F with  F(0)  =1  is  said  to  be  NBU  if 
1 F(t)  for  all  s , t ^ 0 . 

F(s) 

Def Inition : 

The  nonincreasing  stochastic  process  {X(t),  t ^ 0}  is  said  to  be 
NBU  if,  with  probability  1, 

P{T  > s + t I X(u),0  < u < s}  < P{T  > t} 

for  all  s , t , a ^ 0 , -..’here  T^  denotes  the  first  time  the  process 
hits  or  goes  below  a . 

Theorem  3: 

If  the  component  processes  are  Independent  NBU  processes  then 
is  also  NBU. 


Proof : 


Suppose  first  that  $ is  binary  and  let  T denote  the  first  time 
the  process  (jifXCt))  hits  0.  Now  consider 


P{T  > s + t I X^(u),0^u^s,i  = 1,  ...,  n}  = 

E[  (}>  ()( (s  + t))  I X^(u),0^u^s,i  = 1,  ...,  n]  . 


Now  it  follows  from  the  definition  of  an  NBU  process  that  the  conditional 
distribution  of  X^(s  + t)  , given  X^  (u)  , 0 <_  u < s , is  stochastically 
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smaller  than  the  distribution  of  X^(t)  . Hence,  from  Proposition  1 
we  see  that 

E[({)(X(s  +t))  ! X^(u),0<_u<_s,t  = l,  n] 

1 E[4>(X(t))l 
= P{T  > t} 


which  proves  the  result  when  $ is  binary.  As  before  we  can  reduce  the 

j 1 if  i (x)  > a 

nonbinary  case  to  the  above  by  defining  ^ (jc)  = . . B 

^ f 0 if  (jc)  a 

Example  3:  A Semi-Markov  NBU  Process 

If  iX(t),t  ^0}  is  a semi-Markov  process  such  that  X(0)  = m , 
P^^^=l,i>0,  Pqq  = 0 and  the  time  in  state  i is  NBU  then  {X(t)} 
is  NBU.  This  follows  from  the  fact  that  NBU  is  preserved  under  convolution. 


Example  4:  An  NBU  Renewal  Shock  Model 

Suppose  that  X(0)  is  fixed  and  that  at  random  time  points  X(t) 
is  decreased  by  an  i.i.d.  nonnegative  amount.  If  the  random  time  points 
occur  in  accordance  with  a renewal  process  whose  interarrival  distribu- 
tion is  NBU  then  it  is  easy  to  show  that  {X(t)}  is  NBU.  Formally 


I N(t)  ) 

X(t)  = Max  ■ C - y X.,0} 

; i=l  " i 


where  X^  are  i.i.d.,  independent  of  {N(t)}  which  is  a renewal  process 


with  an  NBU  interarrival  distribution. 


Example  5:  An  NBU  Poisson  Shock  Model 


See  Barlow-Proschan  [1],  p.  160. 
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4.  APPLICATIONS  TO  STOCHASTIC  NETWORKS 

The  preceding  theory  can  be  used  in  the  study  of  stochastic  networks. 

For  instance  consider  a network  and  number  the  arcs.  If  we  let  denote 

tlie  capacity  of  the  i—  arc  then  <)>  (x)  could  be  set  equal  to  such  quantities 
of  interest  as 

(i)  4>^(x)  = niax  flow  from  s to  t 

(ii)  l>9(x)  = niax  min  x. 

P ieP  ^ 

wliere  a and  t are  two  given  nodes  of  the  network  and  P denotes  the 
set  of  all  paths  from  s to  t . Thus  equals  the  maximal  flow 

that  can  be  sent  from  s to  t subject  to  the  arc  capacities,  and  <()2(x)  > 

the  maximal  flow  along  a single  path  (previously  studied  by  Barlow  [2]).  If 

we  interpret  x.  as  a cost  of  sending  a unit  flow  along  the  i—  arc  we 

could  be  interested  in 

4>.,(x)  = min  I x. 

P ieP  ; 

the  minimal  cost  path  from  s to  t . Clearly  there  are  many  possible 

functions  41  of  potential  interest,  and  as  it  is  reasonable  to  suppose  ' 

tliat  arc  capacities  (or  costs)  could  change  stochastically  with  time  the  I 

results  of  this  paper  should  be  applicable.  [For  (p^  it  may  be  reasonable  ; 

to  suppose  that  the  x^  increase  over  time]. 

« 

I 


i 
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